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The B~ meson [l| is composed of a bottom quark, b, 
and an anti-charm quark, c, the heaviest quark flavors 
expected to form mesons. The presence of two relatively 
heavy quarks is unique in the system, and affects 
the theoretical calculation of the decay properties and 
mass of the meson. Either quark in the B^r meson 
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can decay weakly, which suggests copious decay modes 
and an expected lifetime much shorter than that of 
other B mesons Q ■ The mass of the meson has been 
predicted using a variety of theoretical techniques. Non- 
relativistic potential models have been used to predict a 
mass of the Bf in the range of 6247 - 6286 MeV/c 2 0], 
and a slightly higher value is found for a perturbative 
QCD calculation [5j. Recent lattice QCD calculations 
provide a Bt- mass prediction of 6304 ± 12 ±^ 8 MeV/c 2 
[6j . Precision measurements of the properties of the B^ 
are needed to test these calculations. 

The Bf meson is too massive to be produced at e + e~ 
colliders operating near the T(4S) center of mass energy, 
and only a few candidate events, which were consistent 
with background processes, were observed at LEP 
With a large 6-quark production cross section and pow- 
erful multipurpose detectors, the Fermilab Tevatron is 
well-suited to the study of all species of B hadrons, in- 
cluding the B^r meson. In particular, this collaboration 
observed the Bf meson in semileptonic decays and used 
these candidates to measure the lifetime of the B^- to be 
0.463^065^0-036 ps *8]. In addition, we have previously 
presented evidence for a Bf — > J/ipn ± signal @ . 

In this Letter, an analysis is presented which observes 
the decay mode B^r — > J/ip n , with J/ip — > /i + /i~ , and 
precisely measures the Bf- mass. This result provides the 
first observation of the B^r meson in a fully reconstructed 
mode. This observation is made in pp collisions at a cen- 
ter of mass energy of 1.96 TeV using the Collider Detec- 
tor at Fermilab (CDF II) [10]. The results presented here 
are based on a data sample with an integrated luminos- 
ity of 2.4 fb , and supercede our earlier measurement 
which was based on an integrated luminosity of 360 
pb" 1 . 

This analysis makes use of the tracking and muon 
identification systems. The tracking system consists of 
double-sided silicon detectors [ll| and a 96 layer open-cell 
drift chamber (COT) [l2j that operate inside a solenoid 
with a 1.4 T field oriented along the beam axis. Charged 
particles originating from the collision point are mea- 
sured in the tracking system with a momentum reso- 
lution of <t{pt)/pt ~ 0.0015p T (GeV/c)- 1 [3. Muon 
candidates from the decay J/ip — * fi + fi~ are identi- 
fied by two sets of drift chambers located radially out- 
side the electromagnetic and hadronic calorimeters. The 
central muon chambers cover the pseudorapidity region 
?y| < 0.6 and are sensitive to muons with transverse mo- 
mentum pt > 1.4 GeV/c. A second muon system covers 
the region 0.6 < |?7| < 1.0 and detects muons having 
Pt > 2.0 GeV/c. Muon triggering and identification are 
based on matching tracks measured in the muon system 
to COT tracks. This analysis is based on events recorded 
with a trigger that is dedicated to the collection of a 
J /tp — > fi + \x~ sample. The first level of the three-level 
trigger system requires two muon candidates with tracks 
in the COT and muon chamber systems that match in 
the transverse view. The second level imposes the re- 
quirement that muon candidates have opposite charge 



and limits the accepted range of opening angle between 
the two muon candidates. The highest level of the J/ip 
trigger reconstructs the muon pair in software, and re- 
quires that the invariant mass of the pair falls within the 
range 2.7-4.0 GeV/c 2 . 

The analysis of the data begins with a selection of well- 
measured J/ip — > fi + fj,~ candidates. Events are required 
to contain two oppositely charged muon candidates that 
satisfy restrictive matching requirements, consistent with 
the full measurement precision of the COT and muon 
chamber systems. We also require that both muon tracks 
have associated measurements in at least three layers of 
the silicon detector and a two-track invariant mass within 
70 MeV/c 2 of the world-average J/ip mass [HI]. This data 
sample provides approximately 17 million events contain- 
ing J /if) candidates, measured with an average mass res- 
olution of 13 MeV/c 2 . 

Both B ± -> J/ipK ± and Bf ->■ J/ip^ combina- 
tions are reconstructed in this analysis. The relatively 
plentiful B ± sample serves as a reference signal and is 
used to develop criteria for the B^r selection . These 
final states are identified by assigning the 7r ± or 
mass to all tracks not used in the J/ip reconstruction. 
In order to provide the measurement resolution sufficient 
to discriminate between directly produced tracks and B- 
hadron decay products, the 7r ± and candidate tracks 
are required to have measurements on at least three lay- 
ers of the silicon detector. Each three-track combination 
must satisfy a fit in which the tracks are required to orig- 
inate from a common vertex and the invariant mass of 
the muon pair is constrained to the world average J/ip 
mass. Approximately 65 000 B^ candidates are identified 
in this loosely selected sample. 

The selection requirements are further improved to 
give a large B^ — ► J/ip signal with very small back- 
ground, in order to increase the significance of the B~- 
observation. Since the B^r lifetime is significantly shorter 
than the lifetime of the B^ , each selection variable under 
consideration has been studied for its effect on preserving 
signal and removing background only for J/ip K combi- 
nations with 80 < ct < 300 /im, where t is the proper de- 
cay time determined from t = Lt • Pt(B) r^ffgpjU , M(B) 
is the mass of the combination, pt(B) is the transverse 
momentum of the combination, and Lt is the transverse 
displacement of the J/ip decay vertex from the beam- 
line. Each selection quantity considered is evaluated for 
its efficiency in retaining the B^ signal and reducing the 
combinatorial background in the 5.4 — 5.5 GeV/c 2 mass 
range. 

Several characteristics of the J/ip candidates are 
considered for selection requirements. Minimum pt re- 
quirements on the and B^ candidates are used to 
suppress backgrounds from J/ip combinations that 
are not related to the decay. Reasonable vertex qual- 
ity is assured by placing a minimum value on the accepted 
probability P(x 2 ) of the mass- and vertex-constrained 
fit used to obtain the B^ candidate. The trajectory of 
the K^- is required to originate from the B^ decay ver- 
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tex by placing a requirement on its distance of closest 
approach dsv(K) an d associated uncertainty <Td sv (K) 
with respect to the vertex found in the J/ip fit. Simi- 
lar quantities dpv(K) and Od PV {K) measured with re- 
spect to the primary vertex are used to remove tracks 
that originate from direct production processes. We sup- 
press the promptly-produced combinatorial background 
by rejecting candidates with small ct. A requirement on 
a c t removes poorly-reconstructed combinations and other 
backgrounds. We also reject J/ip combinations that 
are inconsistent with having originated from the beam- 
line by requiring a small magnitude for the transverse 
impact of the candidate, dpv(B) = hp X pp{B) /\pp{B)\, 
and a small angle j3 between Lp and pp(B). 

Two sets of selection criteria based on these quantities 
are listed in Table |TJ The first step of the selection pro- 
cess is to impose the "standard" selection requirements 
listed above and retain all J / ip K combinations that sat- 
isfy them. However, the selection studies indicate that 
a substantial number of B ± — > J/ipK^ signal events 
fail exactly one of the standard selection requirements. 
To increase the overall signal efficiency, we introduce a 
"high-pr" selection, comprised of combinations that sat- 
isfy more restrictive pp requirements, and pass all but 
one of the other high-px selection criteria in Table [U 
Figure Q] shows the reconstructed B^ — > J/ip sig- 
nal for candidates that satisfy the standard or high-py 
selection criteria. The additional combinations retained 
by the high-pT selection increase the B^ yield by 28%. 
The total signal of approximately 21 100 i? ± candidates 
with a small background of 430 events in the B^ side- 
band region between 5.4 and 5.5 GeV/c 2 demonstrates 
that this selection effectively removes the background to 
B-hadron candidates. 

TABLE I: Selection variables and requirements for the stan- 
dard selection and high-pr selection as described in the text. 
Here "Trfc" refers to the track combined with the J/ip and 
may be a if ± or -n^ candidate for the B^ 1 or Bf respectively. 



Selection variable 


Standard 


High-pT 


p T (Trk) > 1.7 GeV/c 


> 2.5 GeV/c 


Pt{J/iP Trk) 


> 5 GeV/c 


> 6 GeV/c 


P{x z ) 


> 0.1% 


> 1% 


\dsv(Trk)\ 


< 100 pmi 


< 80 pm 


\d PV (Trk)\l '(J dpv (Trk) 


> 2.5 


> 3 


\dpv(B)\/<Td PV (B) 


< 2.5 


< 2 


ct 


> 80 fim 


> 100 /mi 




< 30 /mi 


< 25 /im 





< 0.4 radians < 0.3 radians 



The candidates used for the reconstruction of the Bf 
are obtained by applying the selection criteria in Table U 
to the sample interpreted as J/ip tt^ 1 . Figure shows the 
Bf — > J/tp 7T^~ candidates under the combined standard 
and high-pT selections. The data are also shown for a 



limited mass region suggested by theoretical expectations 

HUH. Aclear 

excess of J / ip 7r ± combinations is evident 
around 6280 MeV/c 2 . 

The data shown in Fig. [DJb) are fitted using an un- 
binned log likelihood function that uses a Gaussian shape 
for the signal plus a background model that includes con- 
tributions from both combinatorial sources and contribu- 
tions from Cabibbo-suppressed Bf — ► J/ip decays in 
the probability distribution function. Simulated events 
are used to estimate the mass distribution expected from 
the misidentified Cabibbo-suppressed process, and it is 
found to be approximately Gaussian, centered 60 MeV/c 2 
below the true B^ mass and having a 30 MeV/c 2 width. 
The characteristic width for the signal in the probabil- 
ity distribution function is varied for each candidate in 
proportion to the mass uncertainty estimated from the 
individual track uncertainties. A similar fit performed on 
the B ± candidates finds that the mass resolution is un- 
derestimated by a factor of 1.55 with this method. Con- 
sequently, the mass uncertainties used in the J/tp ^ fit 
are scaled by this factor. For the central value of the fit, 
the Cabibbo-suppressed contribution is fixed to 0.05 of 
the total Bf yield, which is comparable to the measure- 
ments obtained for Cabibbo-suppressed B^ decays fl4j |. 
Fig. HJb) shows a projection of the fit overlaid on the 
data, which has y 2 / DoF = 34.6/32 and a probability of 
P( X 2 ) = 0.30. 

The result of the unbinned fit over the mass range 
6150 - 6500 MeV/c 2 gives a Bf signal of 108 ± 15 
candidates with a mass of 6275.6 ± 2.9 MeV/c 2 . This 
fit was repeated over the same mass range with the con- 
straint that no signal is present in order to obtain a sig- 
nificance measurement. The ratio of likelihoods between 
the two fits yields a probability for the null hypothesis 
of 1.2 x 10~ 16 , equivalent to a fluctuation of 8a or more 
in a Gaussian distribution. A second significance test 
was also performed, which used a Monte Carlo simula- 
tion to estimate the probability for a background-only 
sample to mimic a signal at least as significant as the 
one observed. This study confirmed the significance es- 
timate obtained from ratio-of- likelihoods test. We con- 
clude that this enhancement is an observation of the pro- 
cess Bf —* J/ipn ± . 

Studies that were performed for the mass measurement 
of other B mesons that contain a J/ip in the final state 
[IH are used for an evaluation of some of the system- 
atic uncertainties. The systematic uncertainty due to 
tracking detector misalignments and material distribu- 
tions used in track fitting is assessed by varying these 
within reasonable limits. Based on these studies, a sys- 
tematic uncertainty of 0.6 MeV/c 2 is assigned to account 
for misalignments and material modeling. A comparison 
between the fitted B ± mass in 12 independent periods 
of data finds the measured mass of the -B meson for 
each subset to be within 1.8<r of the entire sample. Con- 
squently, no systematic uncertainty is assigned for time 
dependent effects on this mass measurement. The mo- 
mentum scale of the tracking system is calibrated with 
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the J/ip, ip(2S), and T states, all of which have been well 
measured previously [3] . The uncertainty of the momen- 
tum scale for the decay products of which was derived 
from this calibration contributes an additional mass un- 
certainty of 0.6 MeV/c 2 . The sensitivity of the mass 
measurement to the assumed Cabibbo-suppressed back- 
ground process was tested with simulated events. Vari- 
ations of 0.0 — 0.1 for the fractional contribution corre- 
spond to mass measurement variations of 0.8 MeV/c 2 , 
which is taken to be the systematic uncertainty associ- 
ated with this unknown background contribution. 

An additional systematic uncertainty comes from the 
fitting procedure. The scale factor associated with the 
mass uncertainty which is used in the fit was set to the 
best value obtained in a fit to the topologically similar 
_B ± — > J/tjj final state. However, other reconstructed 
B hadrons that include a J/ip in the final state have been 
studied as well. Scale factors in the range of 1.25 - 2.0 
have been suggested by these studies, so the fit for the 
B^r mass was repeated for this range. An uncertainty of 
2.2 MeV/c 2 in the calculated Bf mass is indicated due 
to the fitting procedure. The combined systematic un- 
certainty is the sum in quadrature of the uncertainties in 
our calibration, tracking, and fitting procedure. The to- 
tal systematic uncertainty in the B^ mass measurement 
is determined to be 2.5 MeV/c 2 . 



In conclusion, we have observed fully reconstructed 
Bf mesons through the decay Bf — > J/'07r ± . A sig- 
nal of 108 ± 15 candidates is observed with a significance 
greater than an 8a fluctuation of a Gaussian distribu- 
tion. The mass of the B^ meson is measured to be 
6275.6 ± 2.9(stat.) ± 2.5(syst.) MeV/c 2 . 
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FIG. 1: The invariant mass distribution of J/xp K ± combina- 
tions. 
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FIG. 2: (a). The invariant mass distribution of J/ip^ com- 
binations, (b) Identical to (a), but in a narrower mass range 
around the theoretically favored Bf mass. The projection of 
the fit to the data is indicated by the curve overlaid on (b). 



